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We investigate the phase diagram of TE1B4, an Ising magnet on a frustrated Shastry-Sutherland 
lattice by neutron diffraction and magnetization experiments. At low temperature we find Neel 
order at low field, ferrimagnetic order at high field and an intermediate phase with magnetization 
plateaus at fractional values M/M sat = 1/7,1/8,1/9... and spatial stripe structures. Using an 
effective S — 1/2 model and its equivalent two-dimensional (2D) fermion gas we suggest that the 
magnetic properties of TmEk are related to the fractional quantum Hall effect of a 2D electron gas. 

PACS numbers: 73.43.Nq, 71.10Pm, 71.20.Eh, 75.25+z 



The magnetic properties of quantum spins with antifer- 
romagnetic (afm) coupling on frustrating lattices have at- 
tracted widespread interest in recent years due to the dis- 
covery of new types of complex quantum ground states. 
Recent examples include the spin liquid dimer phase ob- 
served in SrCu2(B03)2 [ll (a 2D magnet on the Shastry- 
Sutherland lattice (SSL) QLsupersolid phases on trian- 
gular and kagome lattices Q and the spin ice state in 
pyrochlore 0]. The interest in such systems is in partic- 
ular in magnetization plateaus at fractional values of the 
saturation magnetization as observed in SrCu2(BC>3)2 @| 
and in fractionalized excitations as discussed for spin ice 
@, in particular the possible existence of a magnetic 
monopole Q. This documents that plateaus and frac- 
tional excitations which arc well known in quantum Hall 
physics also emerge in quantum magnets. 

A relation between magnetization plateaus in 2D quan- 
tum spin systems and the quantum Hall effect was estab- 
lished by a mapping of 2D S=l/2 magnets to the 2D elec- 
tron gas Q , requiring an additional Chern-Simons field, 
a fictitious nonlocal strong magnetic field. When this 
field is spatially averaged it combines with lattice induced 
minigaps (known from the 'Hofstadter butterfly' Q) to 
induce fractional magnetization plateaus. In simulations 
using this approach fractional magnetization plateaus 
down to values as small as M sat — 1/9 have recently been 
reported [h| together with complex modulated spatial 
structures of the corresponding ground states. For inves- 
tigations of these phenomena in SrCu2(BC>3)2, however, 
the very high critical field in this material is a severe limit 
and complex spatial structures have not been found ex- 
perimentally so far. 

In the following we present the results of magnetiza- 
tion and neutron diffraction experiments for T1UB4, an- 
other 2D magnet on a SSL and fully accessible to exper- 
iments up to the saturation field. We suggest that this 
compound is an excellent realization of a frustrated mag- 



net on the SSL, we find the emergence of fractionalized 
magnetization plateaus strikingly similar to the case of 
SrCu2(BC>3)2 and we observe spatial structures similar to 
those described in Ref. [nj We find it attractive to spec- 
ulate, on the basis of experimental results to be described 
in the following, that TmB4 represents a magnetic ana- 
logue to the fractional quantum Hall effect (FQHE). 

T111B4 belongs to a group of rare earth tetraborides 
(REB4) that crystallize in a tetragonal lattice (space 
group 127) [Til ] . The rare earth moment position can 
be mapped to a SSL with nearly equivalent bond lengths 
between nearest neighbors. In contrast to SrCu2 (603)2, 
however, (i) T111B4 is a metal and (ii) the Tm 3+ moment 
is J = 6 and subject to strong crystal field (CF) effects. 
Thus the RKKY interaction is expected to be important 
and the basic spin model is an effective S=l/2 model 
close to the Ising limit rather than a Heiscnbcrg model. 

Experimental: Single crystals of TmB4 were grown by 
an inductive, crucible-free zone melting method. The 
residual resistivity ratio was larger than 100, document- 
ing the high sample quality. The neutron sample had 
cylindrical shape with a diameter of ks 7mm and a 
length of re 20mm, isotope enriched n B was used to 
reduce the absorption. For magnetization experiments 
oriented samples were cut with approximate dimension 
lxlxl mm 3 . Magnetization data were measured using a 
commercial magnetometer [l2j ]. For neutron diffraction 
experiments we have used the E10 diffractometer at the 
HMI A wavelength of 0.14 nm was obtained from 

a Ge (220) monochromator. The A/2 contamination was 
removed by Sapphire- and PG-nltcrs. 

Macroscopic properties: T111B4 shows a rich phase di- 
agram as a function of field and temperature (Fig. QJa)) 
14, [lj|. In zero field magnetic order sets in at Tjv = 



11.8K, the afm low temperature Neel phase is stable 
below 9.8 K. Based on magnetization experiments, at 
low temperature three phases are seen: The 'high' field 
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FIG. 1: (a) The phase diagram of TmB4 as derived from 
magnetization data for B || (001) (full symbols), (b) Mag- 
netization at T — 2K (open squares) and T = 3K (open 
circles). The arrows indicate the direction of the field change. 
The inset in (b) gives an overview towards high field, lines 
are measured data, full circles are calculated using the high 
temperature approximation. The magnetic field data are cor- 
rected for demagnetization. 



phase is ferrimagnetic with a plateau at M/M sat = 1/2. 
In the intermediate phase, magnetization plateaus with 
M/M sa t = 1/7, 1/8, 1/9 .. . are seen at temperatures be- 
low 4K. Their observation is subject to hysteresis. After 
zero field cooling plateaus are observed only when the 
ferrimagnetic phase has been reached once, but then in- 
dependent on the direction of the magnetic field changes 
and the field direction. The value of plateau magnetiza- 
tion varies between different runs: although the energy 
difference between plateaus may be very small, the time 
required for a significant spatial rearrangement of mo- 
ments likely exceeds the time available experimentally. 
We note that close to Tjy the intermediate phase splits 
into more complex phases as identified in specific heat 
and resistivity measurements (TBI . Hfl | , but with little sig- 
nature in neutron diffraction and magnetization experi- 
ments. 

The single ion magnetic ground state is derived from 
the magnetic specific heat |16| combined with magneti- 
zation data: The specific heat can be analyzed in terms 
of a CF Hamiltonian H = -gfx B BS z + DSl where D < 
is the CF anisotropy. Towards high temperature, the en- 
tropy reaches the value R In 2 which suggests an effective 
doublet ground state for the noninteracting Tm 3+ ion. 
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FIG. 2: The magnetic structure of Ising spins in TmEU in 
the tetragonal plane. Only domains along the a-direction are 
shown. Red(blue) circles indicate up(down) moments, the 
different colors indicate the periodicity of the various stripe 
structures. Labels (i) to (iii) refer to the phase labels as used 
in Fig. QI a). In (i) the position of the Tm 3+ ions in the 
tetragonal plane is expanded with the SS bonds marked in 
red. Note the almost identical length of exchange paths Ji 
and J2 . 



At low temperature, for magnetic field parallel to the 
c— axis, the magnetization reaches the saturation mo- 
ment of Tm 3+ for B > 5T. Therefore, the local single ion 
ground state is the doublet Mj = ±6, implying that the 
degeneracy of the J = 6 multiplet is completely lifted by 
the CF. The splitting between ground state and first ex- 
cited state can be estimated from the Schottky anomaly 
flil ]. yielding D = 11 A'. Consequently, the first excited 
doublet with Mj — ±5 is at « 100A and strong Ising 
anisotropy is expected. However, transverse exchange 
will have to be added to the dominant Ising coupling to 
improve the agreement of magnetization and specific heat 
data to results for an ideal Ising model. 

Magnetic structure: The low temperature - zero field 
structure is the Neel antiferromagnet (Fig. (Hi)) with 
moments along the (001) direction, implying ferromag- 
netic (fm) correlations along the Shastry-Suthcrland (SS) 
bonds. This follows from the observation of (100) and 
(300) reflections and the absence of af intensity along 
(00L). 

The diffraction data for the high temperature phase 
(Fig. [2In a)) close to Tjy at B = agree well with an 
amplitude modulated structure of fm stripes parallel to 
the a-axis, separated by a tt domain wall. Any 'checker- 
board' arrangements in diagonal blocks can be excluded. 

For neutron investigations of the intermediate phase 
field cooling was used to avoid the hysteresis. This 
leads to the coexistence of magnetic structures with 
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FIG. 3: (a) and (b) show the neutron diffraction pattern in the 
tetragonal plane for the intermediate and the plateau phase, 
respectively. The figures share the intensity color code given 
in (a). For the indexing of the afm peaks see the text. The 
powder lines arise from aluminum in the 'Orange' cryostat (a) 
and the 'VM3' magnet cryostat [TJ used for experiments in 
a field (b). 

different unit cell size: Fig. EJb) shows a diffrac- 
tion pattern with reflections that may be indexed by 
(2m/7,0,0), (2m/9,0,0), (n/7, K,0) . . ., where m,n are 
integers. Lines of weak "harmonic" reflections result from 
stripe structures parallel to the (If 00) or (0K0) direction 
with a size 7 or 9 unit cells (Fig. [^ii b)), directly related 
to the magnitude of the observed magnetization plateaus 
with M/M sa t = 1/7, 1/9. In addition, the afm reflections 
known from the zero field Neel phase have significant in- 
tensity, giving evidence for fm correlations along the SS 
bonds. 

The ferrimagnetic high field phase shows again peaks 
indexed by 1/8, indicating a large unit cell. Using the 
experimental result M/M sat = 1/2, it can be analyzed 
similarly to the intermediate phase. Here the fm cor- 
relations along the SS diagonals arc changed to afm for 
1/2 of the large magnetic unit cell (Fig. EJiii)) whereas 
the remaining moments are just fm. This again results 
in a stripe structure, any diagonal arrangement of the 
moments can be excluded. 

To summarize the experimental results we conclude 
that the spin interactions appear to be (i) antifcrromag- 
nctic, (ii) highly anisotropic and (iii) governed by the 



Ising like CF anisotropy. They lead to stripe structures 
and magnetization plateaus that resemble the predictions 
for SrCu 2 (B0 3 ) 2 . 

Effective S = 1/2 model: For the interpretation of 
the experimental data at low temperatures we use an 
effective S — 1/2 Ising model on the SSL with afm ex- 
change J\ > on the square bonds and J 2 > on di- 
agonal bonds. In the ideal Ising limit (zero transverse 
exchange) the ground state of the model is the Neel state 
for J i > J 2 /2 and an Ising dimer state (opposite spins 
on the SS diagonals) for J\ < J 2 /2. Starting from the 
Neel ground state, an applied field induces a 1/2 magne- 
tization plateau for 2Ji — J 2 /2 < /jbB < 2Ji + J 2 /2 and 
saturation at still higher fields with discontinuous transi- 
tions between these states. From the data on the extent 
of the Neel and 1/2 plateau phases we conclude J 2 w J\ , 
in agreement with the exchange paths expected from the 
geometric structure. An additional transverse exchange 
(exchanging effective spin components +5 and — |) will 
broaden the discontinuous transitions in an applied field. 
Transverse exchange, however, will be small due to the 
strong single ion anisotropy and numerical calculations 
of the magnetization for a 16 spin SSL show that the 
broadening is below the experimentally observable limit 
for an xy-anisotropy e < 0.1. 

For an understanding of the rich superstructures and of 
the fractional plateaus observed in the regime intermedi- 
ate between afm and ferrimagnetic order, the Ising limit 
is not sufficient (the spatial period in this regime is just 
doubled compared to the period of the Neel phase) and 
in fact several higher order corrections may contribute to 
the observed features such as additional exchange inter- 
actions between further neighbors resul ting from RKKY 
interactions and magnetoelastic effects [l8j. 

Whereas RKKY interaction effects appear as the natu- 
ral choice to describe complex spatial structures in a rare 
earth compound and have, in addition to Ising exchange, 
in fact been used to describe e.g. a variety of spatial 
structures in the compound CeSb [ljj[, this will require 
very special interaction ratios over large distances. Con- 
sidering the similarity to SrCu 2 (B03) 2 this nongeneric 
description is not satisfying and we find it worthwhile 
to consider as alternative possibility the influence of the 
transverse exchange usingthe concepts applied to mag- 
nets on the SSL before [8j. These concepts in fact do 
provide a much more natural approach for an under- 
standing of the experimental findings in T111B4. In the 
following we thus speculate that TmB4 could be consid- 
ered as equivalent to a 2D spinlcss fcrmion gas in a very 
strong magnetic field. This will lead to a proposal for the 
understanding of plateaus and the accompanying spatial 
structures in parallel to the corresponding interpretation 
of the high field phenomena in the well investigated com- 
pound SrCu 2 (B0 3 ) 2 0, EE E3|- This interpretation is 
based on (i) the mapping of spins 1/2 to hard core bosons 
and (ii) the subsequent mapping of hard core bosons to 
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fcrmions. In 2D, step (ii) requires the introduction of 
a fictitious magnetic field, the Chern-Simons field i?cs 
with vector potential a(r) and flux $(r) determined by 
the fermion density p(r), given by V x a = <p p(r)e z and 
$(r) = 2np(r). Here, O is the flux quantum and $(r) is 
the flux through a plaquette adjacent to lattice point r 
(choose once and for all one of the 4 plaquettes adjacent 
to a point r with the fermion density p(r) = Sf. + w). 
Apart from this fictitious magnetic field Hcs there is 
an effective fermion hamiltonian with band width deter- 
mined by the transverse exchange and interaction deter- 
mined by the Ising part of the magnetic exchange J z . 

If the nonlocal vector potential attached to each 
fermion is treated in the average flux approximation as 
described in Ref.[§], the problem is reduced to a gas of 
interacting spinless fermions in a strong magnetic field 
characterized by an average flux related to the magneti- 
zation M of the original magnetic model 



$ = 2ir(p(r)) = 2tt (S z ) + 



1 



= 2tt M 



A noninteracting electron gas on the square lattice has 
been studied starting from a tig ht binding band Q as 
well as from Landau levels [2Ct 21, [22], f 



both ap- 



proaches showing the emergence of a fractal structure 
due to the periodic lattice potential, the Hofstadtcr but- 
terfly Starting from Landau levels leads more directly 
to the quantum Hall effect: Landau band splits into p 
subbands, minigaps result which should be observable in 
additional plateaus of the Hall conductance. Whereas 
a direct experimental observation is not possible due to 
the extremely high magnetic fields required, these mini- 
gaps may have been observed experimentally using su- 
perlattice structures [24J. In our case the relevant mag- 
netic field is not the physical field but the fictitious field 
from the Chern-Simons term which produces magnetic 
flux of the order of a flux quantum per site. It has been 
shown in Ref. Q that the conductance plateaus result- 
ing from the combination of discreteness and fictitious 
Chern-Simons field can be identified as some of the mag- 
netization plateaus observed in SrCu2(B03)2- This was 
done by using the SSL instead of the square lattice and 
treating the electron interaction in mean field approxima- 
tion. The present case is similar in that the SSL applies 
also for T1T1B4, but is different in that the electron in- 
teraction is strong due to the dominating Ising term in 
the magnetic hamiltonian. Thus the situation is that of 
electrons on a SSL in the regime of the FQHE and the 
plateaus of interest will result from the combined effect of 
the lattice (leading to minigaps) and of electron-electron 
interactions (leading to the gaps of the FQHE), a com- 
bination which is expected to result in a large variety of 
possible plateaus. We therefore speculate that T111B4 is 
a magnetic analogon of the FQHE on a periodic lattice. 

Our investigations have shown that the Ising like 2D SS 
magnet T111B4 is characterized by unusual magnetization 



plateaus at small fractional values (1/7, 1/9 ... ) , tied to 
complex spatial structures. We have speculated that a 
consistent description of plateaus and spatial structures 
may be obtained from a mapping of this model to the 
FQHE in a 2D electron gas with strong Coulomb in- 
teractions on a lattice in high field. Whereas alterna- 
tive interpretations should not be excluded at the mo- 
ment, our speculation, if true, would open up the new 
field of studying quantum states similar to the FQHE in 
magnetic materials. We hope that this work stimulates 
further investigations although, in view of the efforts to 
relate microscopic numerical calculations for interacting 
electrons to the FQHE j2f| , the numerical work required 
may be beyond the present possibilities. 
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